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ABSTRACT 



Aims. High angular resolution images of the complex gravitational lens system RXS Jl 131-1231 (a quadruply imaged AGN with a bright 
Einstein ring) obtained with the Advanced Camera for Surveys and NICMOS instruments onboard the Hubble Space Telescope are analysed to 
. determine the lens model and to reconstruct the host galaxy. 
Methods. The lens model is constrained by the relative astrometric positions of the lens and point-like images, and by the extended lensed 
\Q , structures. The non-parametric light distribution is recovered in the source plane by means of back ray-tracing. 

■ Results. 1- Precise astrometry and photometry of the four QSO lensed images (A-D) and of the lensing galaxy (G) are obtained. They are 
found in agreement with an independent study presented in a companion paper. The position and colours of the X object seen in projection 
. close to the lens are found to be only compatible with a satellite galaxy associated with the lens. 
I ■ 2- The Singular Isothermal Ellipsoid plus external shear provides a good fit of the astrometry of images A-D. The positions of extended 
O | substructures are also well reproduced. However an octupole (m = 4) must be added to the lens potential in order to reproduce the observed 
+3 lens position, as well as the Ib/Ic point-like image flux ratio. The ellipticity and orientation of the mass quadrupole are found similar to those 
^ ' of the light distribution, fitted by a Sersic profile. The lens (z = 0.295) is found to be a massive elliptical in a rich environment and showing 
. . possible evolution with respect to z = 0. 

^ ■ 3- The host galaxy (z = 0.658) is found to be a substantially magnified (M ~ 9) luminous Seyfert 1 spiral galaxy. The angular resolution is 
J^j sufficient to see regions where stars are intensively forming. Interaction with a closeby companion is also observed. 
, 4- Finally, in the case of RXS J 1 1 3 1- 123 1 , extended lensed structures do not help much in constraining the lens model. 

_ _ 1 
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1. Introduction 

RXS Jl 131-1231 (hereafter Jl 131) was discovered by Sluse 
et al. (2003) as a quasar located at z s = 0.658 and quadruply 
imaged by a bright elliptical galaxy located at zi = 0.295. A 
thorough study of the observed flux ratios at different epochs 
and in several bands has been presented in Sluse et al. (2006; 
paper I). 

The QSO host galaxy is also lensed and appears as a nearly 
complete (~ 305°) bright Einstein ring, clearly seen not only 
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in the near infra-red, but also at visible wavelengths (angular 
radius 6e ~ 1.8")- Thanks to the "proximity" of the source, the 
ring is much brighter than in other known lens systems showing 
extended structures, such as MG0414+0534 (Falco et al. 1997), 
PG1 1 15+080 (Impey et al. 1998), ER0047-2808 (Wayth et al. 
2005, Dye and Warren 2005) or FOR J0332-3557 (Cabanac 
et al. 2005). Impressive images obtained with the Advanced 
Camera for Survey (ACS) and NICMOS onboard the Hubble 
Space Telescope (HST) are presented and analysed in this pa- 
per. Besides the Einstein ring, they reveal detailed arcs, ar- 
clets and substructures in the ring (see Fig0, corresponding 
to patchy emission in the host galaxy. Since Kochanek et al. 
(2001) state that "the shape of an Einstein ring accurately and 
independently determines the shape of the lens potential and 
the shape of the lensed host galaxy", such images should thus 
represent a potential golden mine to constrain the lens model 
at best as well as to recover the morphology of a QSO host at 
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z = 0.66 with an unprecedented angular resolution. Let us first 
investigate in this introduction what we can really expect. 

1.1. Constraining the lens model 

It is well known that constraints provided only by the relative 
astrometry of the lensed images' of a single point-like source 
are not sufficiently numerous to avoid degeneracies in the lens 
potential, which can then lead to wrong estimates of the H 
value derived from time delay measurements (for a general dis- 
cussion on degeneracies, see e.g. Saha 2000, Gorenstein et al. 
1988). Concentrating on the mass distribution, three kinds of 
degeneracy can be identified: the mass-sheet degeneracy (e.g. 
Saha 2000), the radial mass-profile degeneracy (e.g. Refsdal & 
Surdej 1994 in the circular case and Wucknitz 2002 in a more 
general case) and the degeneracy between the internal ellip- 
ticity, the external shear strength and their relative orientation 
(Witt & Mao 1997, Keeton et al. 1997). 

How can extended lensed structures help in avoiding those 
degeneracies? First of all, since the mass-sheet degeneracy im- 
plies a simple rescaling of the lens mass without affecting the 
image structure, the presence of an Einstein ring will not break 
it, as pointed out by Treu and Koopmans (2002) in their analy- 
sis of PG1 1 15+080. Indeed, this degeneracy would be broken 
with lensing data only if sources at different redshifts were si- 
multaneously lensed (because the critical surface mass density, 
and thus the scaling, changes with redshift). 

Second, the mass-profile degeneracy. Knowing the exact 
value of the slope j3 of the mass distribution (a: oc r~^) is 
required to derive the correct value of H . Arbitrarily fixing 
P = 1 (i.e. isothermal distribution) may lead to a systematic 
error of 10% (Wucknitz 2002). Despite the fact that the ob- 
served f3 value may be affected by the mass-sheet degeneracy 
(Wucknitz & Refsdal 2001), it should thus be determined from 
the lensing data. Extended structures may help at this point be- 
cause astrometric constraints alone are unable to do so, even 
when the number of degrees of freedom is sufficient. Indeed 
any astrometry will be identically reproduced with no exter- 
nal shear and a sheet of constant, critical density (J3 = 0). A 
second point-like source at least, or extended structures is/are 
necessary to raise this pathology. Now, in view of testing the- 
oretical mass radial profiles (such as the halo profile predicted 
in the Cold Dark Matter by Navarro, Frenk and White 1996, 
1997) one should bear in mind that gravitational lensing (GL) 
can only constrain the radial profile within a ring defined by 
the distances of the most external and the most internal im- 
ages (Kochanek 2005). In this context, the cusp configuration 
of J 1 1 3 1 is more favourable than the more symmetrical cross 
configuration (like H1413+1 17 for example) but in such a case 
the extended Einstein ring lies in the radial interval already 
probed by the point-like images (see Fig.0 and thus, does not 
help much. The ideal probe would in fact be a second point- 
like source, located just behind the lens, and giving rise to a 
smaller Einstein ring radius because of geometric considera- 

1 Observed flux ratios must be corrected for differential extinction, 
delayed intrinsic variability and/or microlensing before being used to 
constrain a lens model. 



tions. However, stellar dynamics also appears as a useful com- 
plement to GL to put stronger constraints on the mass profile 
(Treu & Koopmans 2004). 

Third, the external shear, the internal ellipticity of the 
mass distribution and their relative orientation can combine 
in several ways to produce a given net potential quadrupole. 
An asymmetric lensing configuration like J 1 1 3 1 (see also 
B1422+231 in Keeton et al. 1997) and additional azimuthal 
constraints provided by extended structures may help in reduc- 
ing this type of degeneracy. 

Finally, two remarks concerning the constraints provided 
by extended structures. First, they are less robust than the as- 
trometric constraints since they rely on the assumed surface 
brightness conservation law. Second, not all the image pixels 
carry independent information 2 . Indeed, obviously, extended 
structures are relevant to constrain the mass model only when 
different images can be mapped onto each other. Thus, a pri- 
ori, the cusp configuration of J 1 1 3 1 is not favourable in this 
context since the three brightest images are merging and can- 
not be compared over a spatial extent much larger than the PSF 
(in contrast with the case of Q0957+561, Keeton et al. 2000). 
In the language of isophotal separatrices (Suyu & Blandford 
2006; see also Kochanek et al. 2001), the strong magnification 
of a very small region in the source plane produces a nearly 
constant surface brightness arc where it is impossible to pre- 
cisely identify isophotes or flux minima. 

Hopefully, in the case of Jl 131, lensed extended safe- 
structures are clearly identified on the HST images and are a 
priori useful as further constraints on the lens model. 

1.2. The lens fundamental plane 

Another important feature of Jl 131 is the fact that the lensing 
galaxy is bright and rather well separated from the lensed struc- 
tures on the HST images. This makes possible a good determi- 
nation of its photometric properties. Furthermore, Sluse et al. 
(2003) have already shown that the lens is an elliptical galaxy 
located at zi = 0.295. Grapping informations on both the mass 
and the light properties of a non-local elliptical galaxy allows to 
probe the evolution of the M/L ratio and the stellar formation 
history through the position of the galaxy with respect to the 
so-called local Fundamental Plane (FP). The latter is a relation 
between the velocity dispersion, the effective radius and the 
surface brightness of elliptical galaxies (Dressier et al. 1987; 
Djorgovski and Davis 1987; J0rgensen, Franx and Kjaergaard 
1995a, 1995b, 1996). 

Isolated field ellipticals, selected on the basis of their mor- 
phology and colours from the Medium Deep Survey up to 
z ~ 0.7, have been compared to the local FP (Treu et al. 2001; 
Treu et al. 2002) and an evolution of the M/L ratio has been 
found in the sense that distant galaxies are more luminous for 
a given effective radius and velocity dispersion. 

More recently, deep spectroscopy with Keck and VLT has 
allowed to probe the evolution of field early type galaxies up 
to z ~ 1 (di Serego Alighieri et al. 2005, van der Wei et al. 

2 Saha and Williams (2001) have shown that the global shape of the 
Einstein ring is encoded into the point-like image time delay ratios. 
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2005, Treu et al. 2005a,b). A common trend is that the evo- 
lution is mass-dependent: more massive galaxies evolve more 
slowly than less massive ones, probably because the latter have 
significantly younger stellar populations. In fact the evolution 
rate of massive field galaxies is similar to that of cluster galax- 
ies (van der Wei et al. 2005) and mass rather than environment 
is found to drive the evolution (Treu et al. 2005a). 

Gravitational lenses offer a completely different selection 
criterium, i.e. based on mass instead of luminosity, and also to 
redshifts up to z — 1 . A first study by Kochanek et al. (2000) 
has been followed by the works by Rusin et al. (2003) and van 
de Ven et al. (2003), which confirm a statistical evolution of 
the M/L ratio in the B and r bands. Since the data presented 
in this paper allow to estimate the velocity dispersion from the 
lens model fitting and the morphological parameters from the 
fit of the observed surface brightness profile, it will be possible 
to check whether the lens of J 1 1 3 1 belongs to the local FP or 
whether significant evolution is present. 

1.3. The host galaxy 

Since only the matching between the lensed images is used to 
constrain the lens model but not the real shape of the source, 
the latter can also be independently reconstructed. The study 
of the host galaxy and its environment at z — 0.658 is thus 
possible. 

Such studies are interesting because host galaxies and their 
environment are closely related to the AGN phenomenon itself. 
The present view of galactic activity requires two components: 
a central engine - a supermassive black hole - and fueling gas 
feeding the engine through accretion processes (e.g. Antonucci 
1993). Even if "black holes have been discovered in every [lo- 
cal] galaxy that contains a bulge^ and that has been observed 
with enough resolution " (Kormendy and Gebhardt 2001), each 
galactic bulge does not harbour an AGN. Therefore, the limit- 
ing condition for nuclear activity seems to be gas fueling. 

From the morphological point of view, the classical picture 
is that the majority of AGN are hosted by early-type, bulge 
dominated galaxies (i.e. 65 - 85% according to Sanchez et al. 
2004, Dunlop et al. 2003) and that disk dominated gas-rich 
galaxies only host small black holes (Kormendy & Gebhardt 
2001), and weak M v > -23 AGN (Dunlop et al. 2003). The 
feeding (through the AGN phase) and growth of the central 
black hole thus seems related to the bulge history rather than 
to the disk gas content. A coherent explanation relies then on 
merging and galaxy interactions (e.g. Di Matteo et al. 2005). 
Besides triggering stellar formation, merging would supply 
fresh gas to the black hole whose activity would be regulated 
by its own energy release and the gravitational potential of the 
spheroid. If such is the case, AGN should statistically correlate 
with the presence of stellar formation in the host and galaxy 
interactions. 

Sanchez et al. (2004) indeed find that roughly 70% of the 
morphologically early-type hosts have blue rest frame colours. 
This is significantly larger than in inactive galaxies at the same 

3 And the black hole mass correlates with the bulge luminosity, not 
with the total host luminosity (e.g. Magorrian et al. 1998). 



redshift (0.5 < z < 1.1). Focusing on the luminosity of the 
[O in] line, Kauffmann et al. (2003) find that highly luminous 
AGN with 0.02 < z < 0.3 are hosted in bulges with young 
stellar populations. However, Grogin et al. (2005) and Sanchez 
et al. (2004) do not observe a significantly higher fraction of 
merging among active galaxies than among inactive ones. The 
latter authors argue that the morphological signature of interac- 
tion may have already disappeared while the bluer colour is still 
being seen. Finally, Dunlop et al. (2003) claim that the apparent 
lack of enhanced merging they observe in active galaxies is not 
a problem since an accretion rate of only 1 M yr~' is enough 
to activate a quasar. Obviously, the debate is still open. 

From an observational point of view, determining the lu- 
minosity profile and colours of the host as well as the signs of 
merging requires high angular resolution. The bright PSF asso- 
ciated with the unresolved AGN makes the task difficult. Even 
with the narrow HST PSF, only relatively low redshift hosts 
have been investigated so far (typically z < 0.5). The (1 + z) A 
cosmological surface brightness dimming joined to the bluer 
rest frame spectra observed in red filters make the contrast be- 
tween QSOs and hosts even worst with increasing z- Thanks 
to NICMOS observations of z 1 QSOs, Kukula et al. (2001) 
claim that hosts statistically follow the same Kormendy rela- 
tion as 3CR galaxies at the same redshift, while at z — 2, only 
the total host luminosity of the host could be determined within 
a factor 2. 

By stretching the host images, gravitational lensing helps 
in revealing host structures located closer to the QSO and does 
improve the resolution in the source plane (although the im- 
provement is limited and not isotropic, see Sect. [3j. The op- 
portunity to recover the host profile at higher redshifts is thus 
offered for the case of J 1 1 3 1 . A larger S/N due to the larger 
number of collected photons is also expected. However, con- 
trary to some believes, lensing does not decrease the contrast 
between the central QSO and the host (the host surface bright- 
ness is preserved while the flux of the point-like images are 
amplified since the latter are not resolved!), so that the noise 
associated with the central PSF residuals is not smaller than for 
unlensed objects. 

A first statistical study of lensed host galaxies with z > 1 .5 
found in the CASTLES gravitational lensing survey has been 
done by Peng et al. (2004). Most of them are luminous and 
some of them also show blue colours possibly betraying star 
formation in the QSO rest frame. 

In the case of Jl 131, located at z = 0.658, the increased an- 
gular resolution offered by the combination of HST+ACS and 
gravitational lensing helps in recovering photometric substruc- 
tures (e.g. star forming regions) in the host as well as possible 
signatures of merging and/or interaction. In order to keep those 
potentialities, the source must be retrieved in a model indepen- 
dent way. This is done with the technique of back ray-tracing. 
This implies no fitting of any source model or source grid in- 
tensity, while saving much computing time. 

In the following, we first describe the observations and their 
reduction in Sec. |3 then we present our technique to recon- 
struct the source in Sect. [3] In Sect. |4] we detail the successive 



4 



Claeskens et al.: RXS J 1 1 3 1- 123 1 : Lens model and source reconstruction 



steps followed to fit the lens system parameters and we give the 
main results. The latter are discussed in terms of the lens and 
host properties in Sect. [5] Finally a summary on the morphol- 
ogy of this complex gravitational lens system is given in Sect. 

El 

Throughout this paper, we use the cosmological parameter 
values Q = 0.3, A,, = 0.7 and H - 50/i5o km/s/Mpc. 

2. Observations, reduction and photometric 
calibrations 

2.1. Observations 

All observations have been obtained with the Hubble Space 
Telescope (HST) in the context of the CfA-Arizona Space 
Telescope Lens Survey (CASTLES). NIR imaging with NIC- 
2 camera has been obtained in filter F160W on 17-Nov-2003 
and ACS images in F814W and F555W have been collected 
with the Wide Field Camera (WFC) on 22 and 24-June-2004, 
respectively. For each filter, several dithered frames have been 
obtained. The observation log is summarized in Tabled 

2.2. Reduction 

The reduction procedure of the NIC-2 data has been described 
Paper I. Two frames have been discarded due to error flags in 
the brightest point-like images, leading to an effective value of 
NExp = 6 in F160W. On the other hand, a systematic error on 
the relative astrometry of 0.003" has been adopted (Impey et 
al. 1998). 

Regarding the ACS images, they have been individually re- 
duced in 3 steps: 

1. Standard reduction (bias and dark subtraction, flatfield 
and gain corrections) has been performed with the STScI 
CALACS pipeline. 

2. Cosmic ray cleaning has been performed on each individ- 
ual frames with an IDL procedure. For each frame a mask 
of bad pixels has been built from cross-correlation with 
the other aligned frames. The affected pixel values have 
been replaced by the median of the four remaining frames. 
Neighbour pixels suffering from charge scattering (i.e. with 
value above the 5 <x noise level) have also had their inten- 
sity replaced by a bilinear interpolation. We checked that 
no pixel has been replaced in the core of the PSFs. 

3. The geometric image distortion of the ACS WFC had fi- 
nally to be corrected because it strongly affects the pho- 
tometry and astrometry of the sources as a function of their 
positions on the CCD. To reach that goal, we applied the 
STSDAS (V3.3) task Pydrizzle (V3.3), which makes use 
of a fourth order polynomial model of the distortion. This 
model is adequate for characterizing the distortion to an ac- 
curacy better than 0. 1 -0.2 pixels over the entire field of view 
(Pavlovsky et al. 2005). However, since J1131 does only 
extend over a few arcsec, we may expect a lower system- 
atic error on the relative astrometry due to the geometric 
image distortion. Indeed, we checked on the DGEO dis- 
tortion residual images corresponding to filters F555W and 




Fig.l. Colour image of RXS Jl 131-1231 resulting from the 
combination of ACS F555W (blue), ACS F814W (green) and 
NICMOS2 F160W (red) images. 



F814W provided by the STScI that, in the restricted detec- 
tor areas where the lens is observed, the dispersion of the 
residual (DX,DY) does not exceed 0.004 pixel around the 
mean value. Given the ACS scale, adopting a systematic er- 
ror of 0.001" on the relative astrometry is thus conservative 
since it is 3.5 larger than the expected error due to distortion 
residuals. 

One should note that interpolation might also produce a 
slight astrometric error but we think that it is included in 
the frame-to-frame measurement dispersion. 

For each frame, a synthetic PSF has been generated with the 
Tiny Tim v 6.3 software (Krist & Hook 2004) at the position of 
the lens centroid. For ACS PSFs, we noticed that Tiny Tim was 
not properly performing the geometric distortion correction, so 
that we processed the distorted PSFs with Pydrizzle, the same 
way as the images. It is worth noting that these final ACS PSFs 
had to be normalized manually. 

A colour picture of RXS Jl 131-1231 made by combining 
the reduced and coadded images in filters F555W, F814W and 
F160W is displayed in Fig.^ Besides the big yellow arc which 
is formed by 3 merging images of the most central part of the 
host, several thin arcs and arclets appear in blue and green 
colours: they are lensed images of extended substructures in 
the host galaxy (see Sect. 14. 5i . 

2.3. Photometric calibrations 

Absolute flux measurements are reported in the Vega magni- 
tude system according to the following relations: 

m ( Fi60W) = -2.5 log[PHOTNU * CR * ^'(Vega)] , (1) 
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Date 


Instrument 


Filter 


NExp 


Exp(s) 


Scale C'/pix) 


RON (e/pix) 


g (e/ADU) 


17-11-2003 


NICMOS-NIC2 


F160W 


5 


640 


0.0759/0.0754 1 


26 


5.4 


17-11-2003 


NICMOS-NIC2 


F160W 


3 


704 


0.0759/0.0754 1 " 


26 


5.4 


22-06-2004 


ACS-WFC 


F555W 


5 


396 


0.050* 


5.3 


2.0 


24-06-2004 


ACS-WFC 


F814W 


5 


396 


0.050* 


5.3 


2.0 



1 NIC2 plate scale along X and Y respectively (see http://www.stsci.edu/hst/nicmos/performance/platescale l 
* the reported ACS scale is obtained after drizzling 



W(F555W,F814W) = -2.5 \og[CR] + m (2) 

where CR is the count rate, photnu = 1.49816e-6 Jy.Sec.DN -1 , 
F y (Vega) = 1043.5 Jy, m = 25.724 for F555W and m = 
25.501 for F814W, as provided on the NICMOS and ACS up- 
dated Web pages. A systematic calibration error is estimated to 
be 0.03 mag for NICMOS from the same Web page, while we 
adopt a formal value of 0.01 mag for ACS. 

Quality flags reveal that the peaks of QSO images A & B 
are saturated on all F555W and F814W images. However the 
saturation is very light because visual inspection reveals that 
the PSF profile is not altered and masking the central pixels 
does not modify the PSF fitting. We are thus confident in our 
flux measurements reported below and made by means of PSF 
fittings. 

3. Image reconstruction in the source plane 

The first source reconstruction algorithms have been developed 
in the radio domain when lensed extended radio structures were 
the only ones to be resolved. LensClean was built on the basis 
of the Clean algorithm by Kochanek and Narayan (1992) and 
has recently been revisited by Wucknitz (2004). On the other 
hand, the LensMem algorithm was based on Maximum Entropy 
reconstruction (Wallington et al. 1996). In these techniques, 
the shape of the source and the lens model are determined in 
two nested fitting procedures. More recently, in their so-called 
semi-linear inversion method, Warren and Dye (2004) replaced 
the source fitting by a matrix inversion. The two former al- 
gorithms proceed in the direction of propagation of light rays 
while the latter goes in the reverse direction (back ray-tracing) 
and, indeed, produces a deconvolved image in the source plane. 
Brewer and Lewis (2005) have followed a more original ap- 
proach based on genetic algorithms. 

Let us assume now that the lens model is frozen. As in 
the algorithm proposed by Schramm & Kayser (1987), our 
method is based on back ray-tracing and does not rely on any 
parametrized shape of the source. The technique of back ray- 
tracing is easy since it simply makes use of the direct lens equa- 
tion: 

to find the (unique) source pixel (x s ,y s ) corresponding to the 
image pixel (x,,y,) given the displacement angle (d x , d y ). The 
source pixel intensity I s is defined as the observed intensity /,-, 
since the surface brightness is assumed to be preserved. Note 



that if several image pixels are mapped onto the same source 
pixel, the median value of the corresponding pixel values is to 
be taken to remove the effect of deviant pixels (e.g. remaining 
after the PSF subtractions). The PSFs located where the unre- 
solved QSO images are formed must indeed be removed since 
they do not satisfy the surface brightness conservation law. This 
source reconstruction algorithm is thus direct and does not need 
any fitting, any regularisation and is non parametric. The draw- 
back is that the convolution process is not included. This is 
not a problem with the present high angular resolution HST 
images, but if the seeing were of the order of the ring width 
(about 1.5" for Jl 131) a preliminary deconvolution of the im- 
ages would then be necessary. One should note however that 
deconvolution usually produces correlation between the inten- 
sities of neighbour pixels and artifacts in the background. The 
latter would partially compromise a clean source reconstruc- 
tion. 

In practice, our source reconstruction algorithm is based 
on the triangular pixel mapping described in Schneider, Ehlers 
& Falco (1992). Each squared image pixel is divided into two 
triangles which are then mapped onto the source plane with 
the lens equation Q and the deflection angle associated with 
the appropriate lens model. For each source pixel, a vector is 
filled with the intensities of all the image pixels mapped onto 
that pixel. We added an original option. It consists in selecting 
image pixels on the the basis of the local magnification value 
M(i, j) in the image plane. Indeed, selecting only image pixels 
with M(i, j) > 1 allows to ignore the flux observed inside the 
critical line, and thus, to reduce the contaminating flux from 
the lens and to ignore the demagnified regions whose mapping 
leads to large noisy pixels in the source plane. This is a key 
point to make the sequential analysis of the lens system pre- 
sented in Sect. 0] Note that M < 1 pixels do not carry crit- 
ical information since their values either are affected by the 
PSF residuals (image A) or correspond to averages over larger 
source regions and do not contain high resolution spatial infor- 
mation on the source. M < 1 pixels are not needed to recon- 
struct the image of the source since each of them has at least 
one corresponding M > 1 pixel. 

The choice of the pixel size in the source plane is not critical 
in the sense that the intensities in the source pixels are not con- 
sidered as parameters (like in other fitting methods), but rather 
as quantities derived from the data (see also Sect. |4}. It is of 
course natural to choose a size smaller in the source plane than 
in the image plane since the whole source has been magnified. 
However, it is useless to consider an extremely small pixel size 
by adopting al/^( scaling for source pixels located close to 
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a caustic. Indeed, first, if the QSO itself is close to the caustic 
(as it is the case for Jl 131), the residuals from the PSF subtrac- 
tions mapped into the source plane will spoil the highly mag- 
nified central parts of the reconstructed host. Second but most 
important, the lensing model is never perfect, so that small but 
resolved photometric lensed structures in the image plane will 
not be mapped exactly at the same place in the source plane. 
However, if the pixel size in the source plane is not too small 
the lensed images will be mapped onto the same source pixel 
and the source substructures will even be enhanced, at the risk 
of being undersampled. It should also be noted that the S/N 
per source pixel is higher when its size is larger, since a larger 
number of image pixels is mapped on that source pixel. We 
made several tests and checked that no gain in resolution was 
obtained with a source pixel scale smaller than half the image 
pixel scale. 

4. Fitting the gravitational lens system 

The lens system J 1 1 3 1 , with its bright lens, its highly structured 
ring and its 4 superimposed PSFs is morphologically complex. 
Simultaneously fitting the PSFs over a bright and variable back- 
ground, the photometric and the deflecting parameters of the 
lens and recovering the real source shape represent an utopic 
task! Instead, a sequential approach has been followed. We 
identified the following steps: 

1. Fitting of a preliminary lens model (e.g. Singular 
Isothermal Sphere (SIS) + shear) by using only the approx- 
imate astrometry of the point-like images (derived from 
gaussian fitting). 

2. Fitting of the PSF positions and intensities. The bright 
background ring is taken into account by simultaneously 
fitting a parametrized source (e.g. a truncated exponential 
disk profile), lensed by the preliminary model. Precise rel- 
ative astrometry and photometry of the point-like images 
may then be derived in each band and PSFs can be removed 
from the original individual frames. 

3. Ring subtraction. Reconstruct the source plane from image 
plane areas located outside the critical line (see details in 
Sect. [3 and stack the resulting frames to get the best S/N 
result in each band. Map the source back in the image plane 
to build the averaged rings and subtract them from the PSF- 
subtracted images obtained in step-2. 

4. Fitting of the lensing galaxy with a convolved photomet- 
ric galactic model (e.g. Sersic profile) on images obtained 
in step-3. Subtract the fitted lens model from the PSF sub- 
tracted images derived in step-2 and stack them to create 
the best image of the lensed host galaxy in each band. 

5. Fitting of the final lens model. At this stage, constraints 
from extended lensed structures can be added to the final 
astrometric constraints. Several lens models can be com- 
pared. 

6. Source reconstruction in the different bands with the best 
lens model. 

We now describe in more details each step defined above. 
Note that all minimizations have been made using the Minuet 
algorithm developed at CERN (James 1994). 



4.1. Fitting a lens model from astrometry 

The astrometric positions of the point-like images yield pow- 
erful and robust constraints on the GL model. Indeed, unlike 
the flux ratios, they are not affected by external or internal per- 
tubations such as differential extinction in the lensing galaxy 
(Nadeau et al. 1991, Jean and Surdej 1998), delayed intrinsic 
variations, microlensing (e.g. Wambsganss 2001) or substruc- 
tures present in the gravitational potential (e.g. Bradac et al. 
2004). Those photometric effects can easily reach tens of per- 
cent while the astrometry of the point-like images is only af- 
fected at the scale of the Einstein ring of the substructures, 
which is of the order of or smaller than the astrometric mea- 
surement errors. 

The weakness of the astrometric constraints is their 
scarcity. After deducting 4 degrees of freedom for the lens and 
source positions, only 6 parameters can be fitted from the astro- 
metric constraints offered by J 1 1 3 1 . This is however more than 
enough to determine the strength (i.e. the Einstein radius 9e) 
and the azimutal anisotropy (either internal [ellipticity e and 
orientation <p e ] or external [shear y and orientation <p y ]) of the 
lens potential. 

In order to avoid problems associated with the computation 
of errors in the source plane close to the caustic, we minimized 
in the image plane the quantity: 



xl = 



xi 



0~x 



y g -yi 



■S 



■v,< - x k \ 2 jy k -yi\ 



CTvk 



CTyJc 



,(4) 



where (x g ,y g ) and (x k , y k ) are the observed positions of the lens- 
ing galaxy and the k th point-like image, with associated errors 
(o- Xg ,o-y g ) and (o- Xk ,o- yk ), and (xi,y,) and (x k ,y k ) are the fitted 
lens and image positions. 

Following Sluse et al (2003), we fitted a simple SIS + y 
model for this preliminary step. 

4.2. Fitting and subtracting the PSFs 

In order to reconstruct the host galaxy in the source plane, the 
PSFs must first be removed. However the PSF flux determi- 
nations are contaminated by the bright Einstein ring (and con- 
versely), especially in the F814W and F160W filters (see dis- 
cussion in paper I). To minimize this problem, the ring is si- 
multaneously modeled as the image of a truncated exponential 
profile, lensed by the SIS + y model. The parameters of the 
source are simultaneously fitted with the PSF intensities and 
positions by minimizing: 



2 _ \ 1 I *>J L, l 



Xe 



(5) 



>J ! 



where /y is the observed intensity in pixel (i,j) with error cry 
and Ijj is the sum of the intensities of the modeled ring and 
of each PSF in pixel (/, j). xl is computed in a small region 
including the PSFs and the inner part of the ring. Note that the 
lens model is not used to determine the PSF positions but only 
the shape of the ring. 
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Table 2. Relative astrometry, flux ratio (with respect to image A) and absolute photometry in filters F555W, F814W and F160W 
of the lensed QSO images A, B, C, D, as derived from PSF fitting (see Sect. 14.21 . After the Einstein ring has been subtracted, 
relative astrometry and absolute photometry are derived for the lensing galaxy G from the best fit of the Sersic model (see 
Sec l4.3> . After subtracting the lens, the astrometry of component X is derived from a gaussian fitting and its magnitude from 
aperture photometry (corrected for infinite aperture). Theoretical magnification ratios Ma/M, and the (unobserved) AGN position 
are derived from our fiducial lens model (see Sect. 14.41 . Magnitudes of the host are calculated in Sect. 14.51 Quoted errors are 
standard errors on the mean, computed from the results on the NExp individual frames, including the systematics estimated in 
Sects. I2~2l and l2. 31 Boldface indicates the final, weighed average astrometry derived from all bands. 
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Table 3. Best values of the Sersic and de Vaucouleurs (dVc) profiles fitted to the lensing galaxy (see Sect. l4.3l for details). Quoted 
errors are the l<x dispersion of the results on the NExp individual frames. The number of d.o.f. is about 3120 for F555W and 
F814W frames, and about 1 100 for F160 frames. 
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The initial conditions are provided by the PSF positions 
given by simple gaussian fittings of the point-like images, by an 
extended circular source and by the SIS + y fitted in Sect. 14.11 
After each individual frame is processed, the astrometry and 
flux ratios relative to image A and the magnitudes of all point- 
like images are averaged on NExp separately in each spectral 
band. The results are reported in Table|2] The quoted errors are 
lcr standard errors (i.e. 68.3%) on the mean. 



The derived astrometric results of the four lensed QSO im- 
ages are mutually consistent between each photometric band. 
This is reassuring and allows us to compute the final weighed 
averaged astrometry also listed in Table [5]in bold characters. 
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4.3. The lens and the X component 

Once the PSFs have been subtracted, most of the light associ- 
ated with the Einstein ring can then be removed as explained in 
step- 3. 

Eventually, only the lensing galaxy G and component X 
are remaining on the frames and the fit of the lens can be 
made without being heavily biased by the ring contribution. Of 
course, it is still necessary to mask out the regions where ring 
residuals are seen, but X is not thought to bias the lens fitting 
since it is very compact, well separated and much fainter than 
the lens. 

We made use of a Sersic profile (Sersic 1968) to describe 
the surface brightness 2 of the lens: 

Ssersic = ^e~ k N ' , (6) 

where r e is the effective radius including half of the total flux, 
S e is the surface brightness at r e and n is the power law index 
(n — 4 corresponds to a pure de Vaucouleurs profile prototyp- 
ing elliptical galaxies, while n — 1 is equivalent to an expo- 
nential profile describing galactic disk profiles) and k is a nor- 
malization coefficient: k = 1.9992n - 0.3271 (e.g. Simard et al. 
2002). Because the isophotes can be elliptical, r is parametrized 
as r 2 (x,y) = x 2 + y 2 /q 2 where q s is the axis ratio of the Sersic 
model. Oxy is the referential attached to the lens principal axes, 
whose positional angle is tp. Before fitting, the Sersic profile 
was convolved with the appropriate Tiny Tim PSF (see Sect. 
I2~2l . 

Minimizing equation [5] with Ijj given by the latter model 
yielded the best parameter values. We also used the Galfit fit- 
ting algorithm (Peng et al. 2002) to check the robustness of our 
results. 

The best model position relative to image A and the best 
model total magnitude are reported for each band in Table |2] 
The best shape parameters are displayed in Table [3] where 
He(z) = -2.5 log Scale 2 is the standard surface magnitude at 
r e {Scale is the pixel scale on the sky) and the physical size of 
the effective radius R e = D&ir e , where Daj is the angular cos- 
mological distance to the lens /. We provide not only the shape 
parameters of the best Sersic model, but also the ones corre- 
sponding to the best de Vaucouleurs model (i.e. when keeping 
n — 4 constant), for further analysis and comparison with pre- 
vious studies of the fundamental plane (see Sect. 15.1.31 . 

The main contribution to the x 2 value (see Table is a 
faint, spatially variable, diffuse background remaining within 
the Einstein ring. This background is very likely associated 
with secondary images of faint, remote, off-axis regions of the 
extended host galaxy. This effect is wavelength dependent since 
the host galaxy is much brighter and extended in the NIR (see 
Sect.l4~5l. 

On the other hand, extended residuals are seen close to 
the lens centroid, in the F814W and F160W bands, at the 4cr 
and 7<t level respectively (for comparison, the X object leaves 
residuals at the 15 and 27 cr level, respectively; see Fig. |2). 
Although the significance is a little bit less in the F814W band, 
the higher angular resolution and the lower background clearly 
reveal a double-object-like signature, more or less oriented 
along the main axis of the fitted Sersic model. 



We made several tests to identify the possible origin of 
those residuals, but without much success. First, they are defi- 
nitely not associated with a putative fifth lensed image, since a 
Sersic-plus-point-like model does not improve the fit. Second, a 
simple bulge+disk decomposition does not solve the problem, 
even whitout co-centering. Third, allowing for boxy isophotes 
(e.g. Peng et al. 2002) does not change the shape of the residu- 
als. Eventually, we simultaneously fit two Sersic profiles re- 
stricted to the most inner part of the lens (0.75" diameter). 
The central residuals are then removed. However, a larger x 2 
is obtained in the wings, the Sersic exponents are very low 
(«i = 1.58 and «2 = 0.82), and most importantly, the results 
cannot be reproduced in the F160W band. Indeed, fixing the 
relative position of the two components to those obtained in 
the F814W band leads to unphysical values of the parameters 
in the F160W band without improving the^ 2 value. 

We also investigated possible artefact origins. For exam- 
ple, the ACS distortion correction and the poorly sampled PSF 
might introduce artefacts during resampling. But the results 
are stable from frame to frame obtained at different dithered 
positions. We recently found out that galaxy-core residuals 
in Sersic-subtracted F814W data have also been noticed by 
Bolton et al. (2006) who get rid of them by applying a more 
complex B-spline galaxy model. 

Finally, we checked that the flux associated with the residu- 
als is smaller than 1% of the total flux of the lens. Lacking high 
angular resolution 2-D spectroscopy, we decided to keep the 
results of the one-component fitting analysis, even if the spa- 
tial extension of the residuals in F814W (roughly along the Oy 
CCD axis) is probably responsible for the slight offset and the 
larger error of the relative position of the lens in right ascen- 
sion, with respect to the other bands (see Table[2}. Differential 
imaging performed after renormalization between the different 
bands confirms the presence of a slight offset of the lens in 
the F814W band. However a map of the F555W-F160W appar- 
ent colour index does not show any fluctuation above the noise 
(0.12 mag rms) which could account for the presence of a dust 
lane. 

A more detailed discussion and analysis of the lens proper- 
ties is postponed to Sect. 15.1.21 

The lens photometric model being subtracted, the last com- 
ponent seen on the residuals is the compact, unresolved X- 
object already identified in Paper I on the NICMOS frames. 
It is also seen on the F555W and F814W ACS images. In each 
band, its position relative to image A has been derived from 
gaussian fits and its magnitude has been obtained from flux 
measurements within a 2 pixel radius aperture corrected for 
infinite aperture. Those positions and magnitudes are listed in 
Table|2 The discussion about the nature of X is delayed to Sect. 

EEU 

4.4. Fitting different lens models 

In order to determine the optimal lens model, we used the best 
available astrometry with a conservative positional error of 3 
mas for each component (see Sect. 14.21 and Table On the 
other hand, making use of the lens photometric model deter- 
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Table 4. Fitted values of the parameters of simple mass models (see text for details). For each model, the first (resp. second) line 
gives the results when fitting (resp. not fitting) the lens position (LF; resp. NLF). 9 E is the angular Einstein radius, s the intrinsic 
ellipticity, y the external shear, 9 E and 9 y the position angles, (A a , A s ) the differences between the modeled and the observed lens 
positions in a and in 5. 
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Table 5. Fitted values of the parameters of best SIE + y model after introducing anra = 4 multipole (£4, 64). The lens position is 
fitted. Same notations as in Table0] 
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Fig. 2. Residuals of the central 2" x 1.5", in F555W (left), 
F814W (center) and F160W (right). The grey scale is in the 
range — 5cr, +5<x (orientation of the F555W images). 



mined in Sect. 14. 31 we built PSF-subtracted and lens-subtracted 
images only displaying the lensed host galaxy. The latter im- 
ages show many extended rafr-structures which may be used as 
new constraints on the lens model (see discussion in Sect.0. 

Those constraints have been combined to the astrometric 
constraints in the minimisation process through the relations 
10} and 0, according to the following relation: 



X 2 = A^/d.o.f a + (1 - A)xt/d.o.f e 



(7) 



where A is a weighing factor, d.o.f. aje are the numbers of de- 
grees of freedom associated with the astrometry and the ex- 
tended image respectively. The value of d.o.f. f is obtained by 
subtracting the number of parameters from the number of pix- 
els, but for the astrometric constraints, we fixed the minimum 
value of d.o.f. fl to 1. Should the number of parameters be too 
large for the number of astrometric constraints, the d.o.f. c is 
more than sufficient. 

For each step of the fit, xl was computed from the recon- 
structed ring intensities Zy, The sum can only be made over pix- 
els (/, j) which are multiply imaged (otherwise Zy - Zy is identi- 
cally zero) 4 . More specifically, the sum was restricted over 68 1 



4 Strictly speaking, the x 2 estimator cannot be applied since the 
computed pixel values are not independent from the observed ones. 
However, the obtained values for the x 2 being not acceptable from a 
pure statistical point of view, a more appropriate statistical test is not 
critical. 



pixels belonging to the sub-structures of the ring well visible 
in the F555W band. 



4.4.1. Canonical lens models 

First, we dealt with three simple models, namely the Singular 
Isothermal Sphere plus external shear (SIS+y, e.g. Schneider 
et al. 1992, Claeskens et al. 2001), the Singular Isothermal 
Ellipsoid (SIE) (Kassiola and Kovner 1993, Kormann et al. 
1994, Keeton and Kochanek 1998) and the SIE+y. We adopted 
the SIE normalization proposed by Kormann et al. (1994), 
where the projected mass within a given surface density con- 
tour is independent of the axis ratio q of the matter distribution. 

The ellipticity e is related to q through q - Jj^ (e.g. Keeton 
and Kochanek 1998). 

Results are given in Table 0] firstly when only fitting the 
astrometry of the 4 point-like images and of the lens (i.e. A — 1 
in Eq.Q LF case). The dominant contribution to thexl comes 
from the lens position, which is not well reproduced by any of 
the models. So, when we added new constraints by including 
extended structures (A = 0.9), we removed the lens position 
from the fit, hence the "no lens fit" (NLF) case in Table|4] 

A striking result is that in the LF cases the astrometric po- 
sitions are better reproduced in term of xl when an external 
shear is present (i.e. with the SIS + y or the SIE + y model). 
However, when the lens position is not fitted (NLF cases), the 
astrometry of the point-like images is much better reproduced 
with the SIE + y model (xl = 10), as generally expected when 
both internal and external sources of anisotropy are included 
in the model (Keeton et al. 1997). However, the offset between 
the predicted and observed lens positions is significant since it 
amounts to about 14 <r. 

The SIE + y model is also slightly favoured by the extended 
structures, in terms of xl- So despite the lens position offset, 
we shall call this latter model the fiducial lens model in the 
following. 
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Let us come back to the constraints offered by the extended 
substructures. The ^/d.o.f. c values listed in Table @] are not 
strongly sensitive to the adopted lens model! Anticipating on 
Fig. [6] where corresponding substructures are correctly iden- 
tified with our fiducial model - meaning a reasonably good 
astrometric matching -, we can check on the direct F555W 
and F814W images that the surface brightness is not identi- 
cal between the corresponding substructures. The hypothesis 
of surface brightness conservation is thus violated . So, even 
a perfect lens model would lead to a significant disagreement 
in terms of xl- The reason why the surface brightness is not 
preserved may be manifold. First, in the case of the present 
HST data of J 1 1 3 1 , the thin arcs are probably not resolved radi- 
ally and could thus be partially amplified. Second, spatially re- 
solved sources still usually consist of many unresolved sources 
on which microlensing might act, resulting in a possible lo- 
cal violation of the surface brightness conservation hypothesis. 
This would especially be true at high angular resolution. Third, 
differential extinction in the lensing galaxy may also affect the 
apparent surface brightness, although, it should not be very ef- 
ficient in the F814W band. All those arguments should con- 
stitute a warning for lens modeling only relying on extended 
structures. 

In this context, the y — e degeneracy is not entirely raised 
by the extended structures. However the asymmetric cusp con- 
figuration of Jl 131 is more constraining than symmetric quads 
(see the case of B 1422+231 in Keeton et al. 1997). 

Finally, let us mention that the source position and theoret- 
ical amplification ratios of the point-like images given by our 
fiducial model are listed in Table |2 As already pointed out by 
Sluse et al. (2003), they do not agree with the observed flux 
ratios. Microlensing is indeed observed in at least one compo- 
nent, but substructures in the lens potential do not seem to be 
favoured (see Paper I for a detailed discussion of this issue). 

4.4.2. Influence of the X-object 

As demonstrated in Sect. 15.1.11 the X-object is very likely as- 
sociated with the lens. Thus, one can wonder whether it may 
influence the deflection potential of the main lens and whether 
it is responsible for the reported offset between the predicted 
mass centroid and the observed light centroid of the lens (see 
Table[3J. Regarding the second question, the X-object can not 
be the cause of the offset, since it is located on the North side 
of the lens (see Fig.[Qand Tabled, while the reported offset is 
essentially in the South direction. 

To answer the first question, we followed two different 
approaches. In the first one, using the lensmodel vl.06 soft- 
ware (Keeton 2001), we modeled the X-object as an SIS, with 
its position fixed to the observed one. Adding this compo- 
nent to the SIS + y model did not improve the xl an d the 
value of the Einstein radius was found to be extremely small 
(#e,x ~ 10~ 5 ")- In the second approach, we fitted several non- 
parametric asymmetric models (searched by a Monte-Carlo 
method), using the PixeLens public software (Saha & Williams 
2004). We failed to predict any mass overdensity in the vicinity 



of X. Thus, both results tend to show that the X-object does not 
play an important role in the lens potential. 

4.4.3. More general lens models 

Although the number of astrometric constraints is limited and 
despite the fact that the extended substructures may not be very 
helpful because the hypothesis of conservation of the surface 
brightness may be violated (see Sect. 14.4. 1> . it is tempting to 
use more complex lens models in order to try to better repro- 
duce the observed lens position. We have tested the Softened 
Power Law Elliptical Mass Distribution (SPEMD) and the ad- 
dition of higher order multipoles to the SIE + y model. 

The deflection angle of the SPEMD can be efficiently com- 
puted with the fastell algorithm provided by Barkana (1998). 
Allowing a variable index of the radial profile and/or a finite 
core radius and proceeding as in Sect. 14.4.11 did not improve 
the fit, which indeed converged toward the singular, isothermal 
profile. 

The lens potential may be written in terms of a multi- 
pole expansion (Kochanek 1991, Schneider et al. 1992). While 
the monopole (m - 0) is related to the radial profile and the 
quadrupole (m = 2) is arising from an elongated mass distribu- 
tion, the octupole term (m = 4) indicates the presence of quad- 
rangularity (e.g. boxiness/diskiness; Trotter et al. 2000), which 
is observed in the isophotes of some elliptical galaxies (e.g. 
Bender et al. 1989, Rest et al. 2001). Such a term in the lens 
potential is also known to influence the magnifications of the 
point-like lensed images and has been introduced in the con- 
text of "anomalous" observed flux ratios (Evans & Witt 2003, 
Keeton et al. 2003). 

Since the observed flux ratio of J 1 1 3 1 are not well repro- 
duced by quadrupolar models (see Paper I and Table 0, we 
had one more reason to incorporate higher order terms in the 
lens potential. The lensmoded vl.06 software (Keeton 2001) 
has been used to simultaenously fit the astrometry of the point- 
like images and of the lens, as well as the observed flux ratio 
Ib/Ic, which is thought not to be affected by microlensing or 
by differential reddening (see Paper I). We found that both is- 
sues of the "anomalous" flux ratio and of the lens astrometry 
could be solved with an octupole (but not with an m = 3 term). 
This result was confirmed by introducing an octupole in our 
modelisation following the definition (with m = 4) : 

<pm = — — cos m{6 - m ) (8) 
m 

and performing a new fitting. The values of the parameters of 
the best SIE + y + octupole model are given in Table [5] Note 
that the amplitude of the perturbation £4 can be related to the 
amplitude 04 of the deviation of the isophotes of the surface 
mass density, assuming that the unperturbed profile is isother- 
mal: 



The discussion of this model is made in the context of Sect. 

EE21 
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Table 6. Best values of parameters of the Sersic profile fitted on the reconstructed host galaxy. 



Filter 


ji e (z) (mag/arcsec 2 ) 


9s 


y(RA.) 


R e (kpc h^) 


n 


^Sersic.tot 


F160W 


18.53 


0.49 


-67.8 


9.8 


1.05 


17.16 


F814W 


21.15 


0.48 


-66.5 


14.4 


1.25 


19.01 




Fig. 3. Colour image of the reconstructed host galaxy of RXS 
Jl 131-1231, combining the ACS F555W (blue), ACS F814W 
(green) and NICMOS2 F160W (red) images. 



4.5. The host galaxy 

The host galaxy is finally reconstructed from the median of the 
stacked PSF-subtrated and lens-subtracted ACS and NICMOS 
images (see Sect. 14.21 and 14. 3> . as described in Sect. [5] The 
adopted lens model is the best SIE + y + octupole described 
in Sect. 14.4.31 The multicolour reconstructed image of the host 
in the source plane is shown in Fig. [3] Let us already note here 
that adopting instead the SIE + y fiducial lens model defined in 
Sect. l4.4.T1 vields qualitatively the same results. 

In order to avoid artifact close to the center, we made aper- 
ture photometry in a ring with inner and outer radii of 0.12" 
and 1.5" respectively. Since the source is reconstructed in each 
photometric band from the unique median image, it is impos- 
sible to derive error bars from the dispersion of the results. 
Rather, we computed the Einstein ring corresponding to the re- 
constructed source (whose regions with r < 0.12" and r > 1.5" 
have been masked) and we subtracted it from each image in the 
image plane. We derived the error on the magnitude from the 
dispersion of the residuals. This technique overestimates the 
noise since many pixels in the source plane have an increased 
S/N thanks to multiple imaging. Photometric results are listed 
in Table|2] as well as the position of the central AGN predicted 
by our fiducial model. The large astrometric error bars reflect 
the uncertainty on the predicted lens position. 

On the other hand, the ratio between the flux of the Einstein 
ring built as described above and the flux of the reconstructed 



source gives the total magnification of the host in each band. 
We found M 555 = 10.9, M 8 i 4 = 9.9 and M m = 7.8. The 
decrease of the magnification from F555W to F160W corre- 
sponds to the fact that the source is more extended in the near 
IR, well beyond the caustic. 

Finally, although the host is quite irregular, we tried to fit 
a Sersic profile to quantify its morphology. We kept the host 
position fixed to the AGN position computed by the lens model 
and we considered a thick ring with 0.12" < r < 1.5". Since 
the noise in the source plane is ill-defined, the^ 2 value and for- 
mal errors are meaningless. We preferentially used the F160W 
image where the host profile is the smoothest, but we checked 
the stability of the results by also fitting the F814W images. 
Results are listed in Table |6] Structural parameters look rather 
stable, although the effective radius and the slope are slightly 
modified in the F814W band. The host galaxy profile is best 
reproduced by a Sersic model with n ~ 1 typical of disk dom- 
inated galaxies (e.g. de Jong et al. 2004). Although the host is 
globally well subtracted, positive residuals are correlated with 
emission in the F555W filter, while negative residuals are seen 
around the center and might indicate that the PSFs have been 
slightly oversubtracted. Total magnitudes of the fitted model in 
a box of 4" x 4" are also given in Table|6] Again, analysing the 
host galaxy reconstructed with the fiducial model yields results 
very close to those presented in Table[6] 

The physical properties of the host are discussed in Sect. 

I5~2l 
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5. Results and discussion 

5.1. The X-object and the lens 
5.1.1. The X-object 

First, although X does not seem to be resolved, it cannot be 
the fifth lensed image of the central AGN. Indeed, besides the 
fact that its colours are completely different from the ones of 
images A-D (see Tabl^5]and Fig.^, a non-singular lens model 
would only predict a faint lensed image roughly aligned with 
G and D. 

Second, X is not a secondary image of a substructure in the 
source. If such was the case, the primary image should have 
been detected since it would be well separated from the ring 
and at least as bright as X, depending on whether it is resolved 
or not. Fig.@]shows that it is not the case. 

Third, X is not a foreground star. Indeed, Table |5] shows 
that its apparent colours are (F555-F814) x = 2.01 ± 0.31 and 
(F814-F160) x = 2.14 ± 0.08. Stars with V - I * 2 are either 
giant K5 or MO dwarfs (Drilling and Landolt 1999). But these 
stars have V-H smaller than 3.6 (Binney and Merrifield 1998) 
and are consequently bluer than the V-H colour index observed 
forX. 

In conclusion, X is very likely a satellite dwarf elliptical 
galaxy of the lens. Indeed, its apparent colours are very similar 
to that of the lens ((F555-F814) c = 1.97 ± 0.14 and (F814- 
F160) G = 1.86 ± 0.10; see Fig.|TJ. Its restframe colours ((B - 
r), = 1.31 ± 0.31 and (B - J) z = 3.50 ± 0.31) are comparable 
to those of local ellipticals (see Sect. 15.1.31 . 



the octupole simultaneously helps in reproducing the observed 
flux ratio I-rHc- 

However the physical origin of this octupole is difficult to 
identify. Indeed, to improve the fits, the octupole position an- 
gle must be completely different from that of the quadrupole. 
This means that the solution (£4, #4) is degenerated with the 
solution (-£4, 64 + 45°) and that the boxiness/diskiness cannot 
be assessed. Using Eq. [9] the surface mass density isophotes 
are characterized by |«4| = 0.06. This value is quite high when 
compared to the boxiness/diskiness of light isophotes observed 
in elliptical galaxies (Rest et al. 2001). Besides that, as we saw 
in Sect. 14.31 allowing for the boxiness of the isophotes did not 
remove the residuals observed in the F814W band. And why 
are those residuals mainly seen in the red bands? The mass 
octupole might thus be a (non-unique) mathematical way of 
accounting for some pertubations in the potential, for exam- 
ple due to interactions with a satellite galaxy (X-object?). This 
puzzling situation could be investigated with 2-D spectroscopy 
obtained at high angular resolution. 

In conclusion, while low-mass halo substructures do prob- 
ably not provide the correct explanation for the observed 
'anomalous' flux ratios of Jl 131 (see paper I), the m = 4 mul- 
tipole can reproduce the observed flux ratios and lens position 
but its physical nature is difficult to identify, as it is for other 
lenses (Kochanek & Dalai 2004). 

Finally, note that the value of the external shear (see Table 
[3l is also found to be quite high. However its direction points 
towards the position (affected by large error bars) of a group 
of galaxies identified at the lens redshift from a red sequence 
analysis by Williams et al. (2006). 



5.1.2. Comparing light and mass distributions 

As already noted by several authors, the pure 1/r 4 de 
Vaucouleurs law does not correspond to the best fitting of the 
surface brightness profile of ellipticals. Although the values we 
obtain for the slope parameter are significantly smaller than the 
canonical n — 4 value, they remain typical of early-type galax- 
ies, which exhibit quite a broad range of values (e.g. de Jong et 
al. 2004, di Serego Alighieri et al. 2005). 

The structural parameters of the best Sersic model of the 
lens are reported for each band in Table [3] Its shape is pretty 
circular, especially in the F555W filter, where the derived ori- 
entation is meaningless. In the red bands, the averaged axis ra- 
tio is q s = 0.936 ± 0.004 and is close to the axis ratio of the 
mass distribution, q m = 0.84 - 0.91, as derived from the fidu- 
cial and best models in Tabs. @] and [5] Similarly, the averaged 
orientation of the light distribution (P.A = -57.8+2.2 degree) is 
comparable to the orientation of the mass model (-50. 1 < P.A. 
< -60.8 degrees). This seems to indicate a good alignment be- 
tween the projected luminous and dark matter distributions. 

On the other hand, we noticed in Sect. 14.4.11 that the mass 
centroid of the fiducial model had to be offset by about 0.04" 
with respect to the photometric centroid to reproduce the image 
positions, unless an octupole is added to the lens model (Sect. 
14.4.31 . The presence of an octupole might be related to the fact 
that photometric residuals are also seen in the F814W band 
after subtracting the best Sersic profile. It is also striking that 



5.1 .3. Lens colours and evolution 

Despite the uncertainties relative to the lens model described 
in the previous section, the monopole term, i.e. the value of the 
Einstein ring radius, remains quite robust. This makes possible 
a deeper investigation of the lens. 

In order to seek for galaxy evolution, the lens must be 
located with respect to the local FR defined for example by 
the Coma cluster FR Following the analysis by van de Ven et 
al. (2003) based on previous works by J0rgensen, Franx and 
KjEergaard (1995a, 1995b, 1996), the Coma FP relation can be 
written as: 

\ogR e = alogcr c + Bfi e + y, (10) 

whereor = 1.24 ±0.07,/? = 0.328 ±0.001 andy = -9.12±0.05 
in the Gunn r filter. In Eq. fllOi . R e is the effective radius in 
kpc, o~ c is the central velocity dispersion in km/s and fi e is the 
surface brightness at r e in mag arcsec~ 2 . 

In a given band, the evolution of the M/L ratio with re- 
spect to local galaxies can be quantified from the FP analysis 
using the following relation (Treu and Koopmans 2004 and ref- 
erences therein): 

Alog(M/L) = -0.4^, (11) 

where Ay = y D b s - y . This relation is independent of the pos- 
sible evolution of the slopes a and B (Treu et al. 2005b). 
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Since one easily shows that overestimating /j e or cr c in Eq. 
^| would hide the lens evolution, we first discuss the determi- 
nation of those parameters. 

i) fi e and the lens colours 

Since the B, (Gunn-)r and J bands redshifted at z = 0.295 
nearly correspond to the observed ones F555W, F814W and 
F160, respectively, it is easy to derive the (B - f) z and (B - J) z 
rest frame lens colour indices. It is also possible to estimate the 
rest frame surface brightness in the r-band // e>r (z = 0) needed 
to locate the lens with respect to the local FP. In the latter case, 
we used the following colour transformation: 

/V(Z = 0) = He,FiUw{Z) + (r z,AB ~ -F814W A b) - (fir - C F8 i 4 w) 

-7.51og(l+z), (12) 

where the first term of the right member is the observed sur- 
face brightness listed in Table |3j the second one is the colour 
term between the redshifted Gunn r filter and the F8 14W pass- 
band, expressed in the AB magnitude system (Oke 1974) to 
measure absolute flux densities; the third term relates Vega to 
AB magnitudes; the last term corrects for the (1 + z) 3 cosmo- 
logical dimming for flux densities expressed in frequency units 
as in the AB system. The colour term has been derived by in- 
tegrating the properly redshifted spectroscopic template of el- 
liptical galaxies (Kinney et al. 1996, Mannucci et al. 2001). c, 
and Cp8i4w have been obtained from Frei and Gunn (1994) and 
the ACS/WFC zeropoints 5 respectively. Similar relations are 
derived to compute the B z and J z rest frame magnitudes. 

We found the following colour indices: (B - r) z = 1.27 ± 
0.13 and (B - J) z = 3.17 ± 0.14. These values are similar to 
the colours of local elliptical galaxies ((B — r) = 1.17, Frei & 
Gunn 1994; (B - r) = 1.2, van de Ven et al. 2003 and (B - 
J) = 3.34 ± 0.12, Mannucci et al .2001; (B - J) = 3.25, Fioc 
and Volmerange 1997), with no sign of significant reddening. 
Adding the fact that the surface brightness is estimated in the 
r band, we conclude that the value // e , r (z = 0) = 21.11 ±0.1 
should not be overestimated by dust extinction. 

ii) o~ c : the Faber-Jackson relation and the lens environ- 
ment 

Lacking a high resolution, high S/N spectrum of the lens, 
the central velocity dispersion <x c , needed in Eq. ^3 must be 
derived either from the Faber-Jackson (FJ) relation (Faber & 
Jackson 1976) or from the velocity dispersion of the best fit- 
ting isothermal ellipsoid, <tsie. Ideally, both estimates should 
be equal. 

Adopting the expression of the FJ relation given by Forbes 
& Ponman (1999), and an intrinsic scattering of Alogcr t . = 
0.05, we found a- c FJ = 251 ±29 km/s. 

On the other hand, cr SIE can be derived from the fitted value 
of the Einstein radius 6e according to the following relation: 

cr sm = 186.4 j W^os km/s, (13) 
3 STScI Web site: http://www.stsci.edu/hst/acs/analysis/zeropoints 



Table 7. Observed evolution of the lens, as a function of the 
underlying assumption on <x c : FJ or GL model with possible 
correction for central velocity dispersion (fsm) an d typical en- 
vironment (Env). 



Hypothesis 


cr c (km/s) 


A log(M/L) r 


GL 


355 ± 


35 


0.01 i 


= 0.06 


GL+fsiE 


309 ± 


28 


-0.1 : 


t0.08 


GL+fsi E +Env 


291 ± 


26 


-0.14 


±0.08 


FJ 


251 ± 


29 


-0.24 


±0.08 



where Dos (resp. Dls) is the angular distance from the observer 
(resp. the lens) to the source. From Sect. 14.41 a robust estimate 
of 6e — 1.84 ±0.01" is obtained and a value ofcrsjg = 355 ± 1 
km/s is derived. 

The normalization factor /si E = cr e /crs IE can be estimated 
theoretically under several hypotheses (e.g. assuming that the 
true profile is really isothermal). Kochanek (1994, 1996) and 
van de Ven et al. (2003) found that the expected value of the 
normalisation factor is compatible with 1, with a typical uncer- 
tainty of 10%. However, recent spectroscopic observations of 
gravitational lenses show that / might rather be smaller than 
1. Indeed, with a sample of five lenses, Treu and Koopmans 
(2004) found / SIE = 0.87 ± 0.08 (RMS) and in the case of 
B2045+265, Hamana et al. (2005) even found / = 0.53. 
Adopting the value / = 0.87±0.08, we derived <r CtGL = 309±28 
km/s. Although the values of cr c CL and <x c FJ are quite different, 
the 1-cr error bars nearly overlap. 

However, cr SIE might also have been overestimated due to 
the presence of external convergence originating in the lens 
environment. By means of simulations, Keeton and Zabludof 
(2004) have shown that typical lens environment can produce 
a fractional overestimate of <xsi E by about 6%. Correcting also 
for this effect would lead to the value cr c CL — 291 ± 28 km/s 
in better agreement with the FJ relation. A fractional increase 
of <xsi E by 23% due to the environment would be necessary to 
reproduce the value of <x c FJ . This would require a dense mat- 
ter sheet with k ~ 0.34, which could be provided by associ- 
ated and/or foreground galaxy clusters. Indeed, Williams et al. 
(2006) have identified two Red Sequences in the field, corre- 
sponding to galaxy groups or clusters in the vicinity of Jl 131, 
one at the lens redshift (offset ~ 37") and one at z ~ 0.1 (off- 
set ~ 137", also seen in X-rays, Chandra OBSID 4814, PI: 
Claeskens, public data). 

Hi) The lens evolution 

We have adopted the structural parameters of the best de 
Vaucouleurs model, i.e. R e = 9.3 ± 0.2 kpc (see Table [5}, the 
derived value of yu e f (z = 0) = 21.11 ±0.1 (see above) and 
several estimates of <x e in the range 25 1 - 355 km/s (depending 
on whether the FJ relation or the GL model is used; see above). 
We then made use of Eqs. (II Oi and (11 U . to locate the lens with 
respect to the local FP and to estimate its evolution in the Gunn- 
r band (see left and right panel of Fig.|5] respectively). 

More quantitatively, the evolution in the Gunn-r band of the 
M/L ratio is reported in Tab.0 The evolution obtained with the 
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Fig. 5. Left: Local fundamental plane in the rest-frame Gunn r-band, defined by the linear fit (J0rgensen et al. 1995a,b; thick line) 
of the Coma cluster galaxies at z — 0.023 (small dots). The lens RXS Jl 131-1231 (triangle) is added to the lenses analysed by 
van de Ven et al. (2003). Right: M/L evolution in the Gunn r-band predicted for a single burst of stellar population occuring at 
z* = 4.3 (mean cluster formation; thick line) down to z* =2.0 (maximum progenitor bias correction;shaded region) or at z* = 1 .9 
(best fit model for the lens galaxies analysed by van de Ven et al. 2003; dashed line). RXS Jl 131-1231 is added as a triangle. 
The error bars on the abscisses are included in the symbol size; the errors on the ordinates reflect the dependence on the assumed 
value of cr c (see Tab. [7}. Figures reproduced from van de Ven et al. (2003). 



normalized and environment corrected lens model best agrees 
with the value A log(M/L) r = -0.14 + 0.03, derived from the 
evolution rate Alog(M/L) r /Az = -0.47 + 0.11 of lenses found 
by van de Ven et al. (2003). 

However, except under the hypothesis that the correct value 
of cr c is derived from the FJ relation, for which a stronger evolu- 
tion is found, no estimate can rule out the absence of evolution 
to better than ~ 1 .5cr level. 

In any case, the apparent evolution of the lens with respect 
to local elliptical galaxies is weak since: 

1 . the rest frame colours of the lens, (B - r) z and (B - J) z are 
similar to the ones of local ellipticals (see above); 

2. the total M/L B within the effective radius is not signifi- 
cantly smaller than in local galaxies. Indeed, first assum- 
ing an isothermal mass distribution with asm = 335km/s 
(after a 6% reduction due to the lens environment, Keeton 
& Zabludof 2004, see above), the total mass inside the ef- 
fective radius is found to be M = 3.9 ± 0.7 10 11 M . 
Second the light: integrating the Sersic model within R e in 
the reference frame B, band and adopting M B0 = 5.47 
(Allen 1999), we find L B = 4.96 + 0.5 lO 10 ^ 2 L and 
M/L B = 7.8+1 h 5 Q (M/L B ) . This is quite close to the upper 
limit observed within the effective radius of local elliptical 
galaxies (3.5 < M/L B h 50 < 8.1, Gerhard et al. 2001). 



In summary, although, within the uncertainties on cr c , the 
FP analysis shows a possible weak evolution of the lensing 
galaxy, the latter looks quite similar to local elliptical galaxies, 
in terms of rest frame colours and M/L ratio. This is not unex- 
pected, firstly because a significant scatter is observed between 
individual objects (e.g. van de Ven et al. (2003) in their study of 
a sample of lensing galaxies). Secondly, massive galaxies like 
this lens are observed to evolve more slowly with redshift (e.g. 
van der Wei et al. 2005). On the other hand, the rich lens en- 
vironment probably helps to reconcile the value of <t c derived 
from o"sie and the one expected from the FJ relation. 

A spectroscopic measurement of the lens velocity disper- 
sion will of course provide a more definitive answer to this is- 
sue. 

5.2. The host galaxy 

Visual inspection of the multicolour reconstucted host galaxy 
(see Fig.0 immediately reveals several striking details: 

1 . Blue spots are easily located. These are associated with the 
blue arcs and arclets seen in the image plane (see Fig.Q. 
The associations are indicated in Fig. [6] We identify one 
quadruplet of arclets, two triplets (with one pair of merging 
images in each of them) and three pairs (for which the iden- 
tifications are made from a circular sub-source extracted 
from the real source). 
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Fig. 6. Associations of the blue arcs and arclets in the image plane with the individual blue emitting regions in the source plane, 
based upon the fiducial lens model (illustrated in F814W). 



2. The blue spots are spread along a spiral arm. Keeping in 
mind that the F555W filter roughly corresponds to the U 
band at z = 0.66, we may infer that they are very likely 
associated with intense stellar formation in the host. 

3. On the left side of the host, a faint, blueish companion is 
clearly seen and is likely in interaction with the host, which 
indeed shows disturbances rather than a smooth profile. 

4. A strong radial colour gradient is seen across the host 
galaxy. It is not due to artifacts related to the PSF subtrac- 
tions, since it is already visible in the Einstein ring of Fig. 
^ The relatively bluer (coded as yellow) central emission 
could be due to the presence of many young stars associ- 
ated with a starburst in a nuclear ring (Knapen 2005). This 
could be tested with 2-D spectroscopy. 

The following results can be deduced. First of all, statement 
number 1 underlines the coherence of our fiducial GL model 
since a correct association between the various arcs and arclets 
is possible. 

Secondly, the simultaneous observation of the AGN phe- 
nomenon, of stellar formation and of galaxy interaction il- 
lustrates the current AGN - black hole paradigm, according 
to which the growth of a black hole through the AGN phase 
occurs when galaxy interactions provide gas inflow towards 
the center (e.g. Hernquist 1989, Di Matteo et al. 2005). Since 
the interactions induce gravitational perturbations, simultane- 
ous starbursts are also expected. However, note that even if in 
the past, companions were found more frequently around active 
galaxies (e.g. Dahari 1984), it does not seem to be systematic, 
nor even statistical (e.g. Sanchez et al. 2004, Knapen 2005 and 
references therein). 

Finally and more quantitatively, the source of the lens sys- 
tem J 1 1 3 1 must be a luminous Seyfert 1 spiral galaxy. This 
result relies on the following line of arguments: 

1 . The host is a spiral, disk dominated galaxy, as derived from 
Fig. 13 and the reasonable fit by a Sersic profile (see Table 



[6] and Sect. 14. 51 . This is corroborated by the rest frame in- 
dex (U-V)j. = 0.78 ± 0.65, derived from a transformation 
of the F555W and F814W measurements similar to Eq.1121 
Such a colour index is typical of local spirals and excludes 
elliptical galaxies at l<x (Fioc & Rocca-Volmerange 1997). 

2. The rest frame absolute magnitude of the AGN is 
My,z(AGN) > -22.9 + 5 log /!5o. This is a robust upper limit 
estimated by adopting the flux of the brightest point-like 
image (i.e. image B) and assuming a total magnification of 
10 (i.e. a conservative value given the macro lens model of 
the cusp configuration). Adopting the Einstein - de Sitter 
cosmological model would even lower the upper limit by 
about 0.5 mag. 

3. One "coherence argument" is that even if only ^ 20% of 
Seyferts are late-type galaxies (e.g. Sanchez et al. 2004), 
disk dominated galaxies are not observed to host AGN 
more luminous than M v ^ -23 (Dunlop et al. 2003). 

4. The rest frame absolute magnitude of the host is found to 
be M y> .(host)= -23.5 ± 0.2 + 5log/i 50 . This is about 1 
magnitude more luminous than the most luminous Seyfert 
galaxy in Smith et al. (1986). However, a significant part 
(i.e. ~ 0.5 mag) of this luminosity excess could be due to 
nebular emission associated with intense star formation ac- 
tivity. Indeed, [Om] ,15007, [Om] ,14959 and H^ fall in the 
F814W filter and could also be responsible for the colour 
gradient mentioned above. After correction for this effect, 
this host galaxy would remain among the most luminous 
Seyfert galaxies. 

6. Conclusions 

This paper was devoted to a thorough analysis of high angular 
resolution HST optical and NIR direct imagery of the complex 
gravitational lens system RXS Jl 1 3 1- 123 1 . 

We have provided precise astrometry and photometry in the 
F555W, F814W and F160W passbands for each component of 
the system, namely the 4 point-like lensed images of the cen- 
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tral AGN, the lensed host galaxy, the lensing galaxy and the 
so-called X-object projected close to the lens. We have also 
reconstructed the image of the host in the source plane in a 
non-parametric way, leading to the view of an active galaxy at 
Z = 0.658 with unprecedented details. 

The summary of our results is the following: 

1. A fiducial lens model consisting of a Singular Isothermal 
Ellipsoid with external shear perfectly reproduces the rela- 
tive astrometry of the point-like lensed images. The shear 
direction is found to be compatible (within large error bars) 
with the position of a galaxy group identified at the lens 
redshift by Williams et al. (2006). However, a disagree- 
ment of ~ 0.04" (i.e. ~ 14cr) is found between the ob- 
served and predicted astrometric positions of the lens. This 
issue could be solved by adding to the fiducial model an 
m = 4 multipole leading to boxy/disky surface mass den- 
sity isophotes. The latter model simultaneously reproduces 
the observed 'anomalous' flux ratio Ib/Ic, confirming that 
lensing by substructures in the halo may not be required 
(see also Paper I). However, the physical nature of the fit- 
ted octupole is difficult to interpret and, as usual, no claim 
is made about the unicity of the proposed GL model. 

2. Faint but significant residuals are found when modelling 
the lens photometric profile in the F814W band with a 
Sersic profile, even when allowing for boxy isophotes. A 
slight offset with respect to the positions in the other pass- 
bands is also observed, although no trace of dust is seen in 
the {B - K) z index. Would this be a trace of past interac- 
tion, maybe due to the X-object? Are both issues 1- and 2- 
related? High angular resolution 2D spectroscopy is proba- 
bly required to answer those questions. It is however strik- 
ing that this kind of residuals in the ACS F8 14W band have 
also been reported for other lenses (Bolton et al. 2006). 

3. The different values of cr c derived from lensing and from 
the FJ relation indicate that the contribution of the lens 
environment could be important. The latter looks indeed 
quite rich (Williams et al. 2006, Chandra observations, PI: 
Claeskens, public data). 

4. Although, within the uncertainties on <x c , the FP analysis 
shows a possible weak evolution of the lensing galaxy, the 
latter looks quite similar to local elliptical galaxies, in terms 
of rest frame colours and M/L ratio. Presence of dust is not 
detected. 

5. The X-object is identified as possibly being a satellite 
galaxy of the lens, but with very little effect on the lens 
model. In any case, it is not the fifth lensed image of the 
AGN nor a lensed substructure of the host. 

6. The reconstructed host is identified as a luminous Seyfert 
1 spiral galaxy showing patchy emission in the UV rest 
frame, probably corresponding to intense star formation. 
A blueish closeby companion is also detected. The global 
shape and magnitudes of the host do not strongly depend 
on the addition of the octupole in the lens potential. 

7. This patchy UV emission is lensed into multiple thin arcs 
and arclets in the F555W band, whose positions are rather 
well reproduced by our lens models. 



8. In the peculiar case of the present cuspy lens system, ex- 
tended structures did not turn out to be very helpful to con- 
strain more general lens models. This underlines the need 
for robust astrometric constraints and this indicates that 
care must be taken when only relying on extended images 
to constrain a lens model. 
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